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We show the existence of a new gapless collective excitation in a rotating fluid system with chiral 
fermions, named as the Chiral Vortical Wave (CVW). The CVW has its microscopic origin at 
the quantum anomaly and macroscopically arises from interplay between vector and axial charge 
fluctuations induced by vortical effects. The wave equation is obtained both from hydrodynamic 
current equations and from chiral kinetic theory and its solutions show nontrivial CVW-induced 
charge transport from different initial conditions. Using the rotating quark-gluon plasma in heavy 
ion collisions as a concrete example, we show the formation of induced flavor quadrupole in QGP 
and estimate the elliptic flow splitting effect for A baryons that may be experimentally measured. 
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Introduction. — Anomalous transport effects in many 
body systems with chiral fermions have generated great 
interests recently. Such phenomena span a wide range of 
physical systems mm, from semimetals to cold atomic 
systems, and from hot quark-gluon plasma (QGP) cre¬ 
ated in heavy ion collisions to cold dense matter in neu¬ 
tron stars. These systems provide possible environments 
with nonzero macroscopic chirality and can manifest the 
microscopic chiral anomaly in macroscopic transport pro¬ 
cesses that would be normally forbidden by symmetries 
like parity invariance. 

One way to induce the anomalous transport effects 
in such chiral system is to apply external electromag¬ 
netic fields. A famous example is the so-called Ghiral 
Magnetic Effect (GME) in which an electric current can 
be generated in parallel to an external magnetic field . 
The GME could lead to experimentally measurable ef¬ 
fects both for the QGP in heavy ion collisions ims] and 
for certain Dirac and Weyl semimetals [16]. Other inter¬ 
esting examples include e.g. the Ghiral Separation Effect 
(GSE) [T7| [TH] , and the Ghiral Electric Separation Effect 
(GESE) For reviews see e.g. 

The anomalous transport effects can also occur when 
the fluid undergoes a global rotation quantified by a vor- 
ticity uj = X V where v is the flow velocity field. 
Such vortical effects were suggested in |3| and quantified 
in holographic models [Ml - ES] and the anomalous hydro- 
dynamic framework m- The so-called Ghiral Vortical 
Effect (GVE) quantifies the generated vector current Jy 
as [21H27| 

Jv = , (1) 

TT^ 

and the generated axial current Ja as pifET] 


where T is temperature and /r and are vector and 
axial chemical potentials. It was suggested that such 
CVE may lead to baryon charge separation in heavy ion 
collisions ^511^ . 

The vorticity-driven anomalous transport effects in 
Eqs. Q and ([^ couple together the vector and axial 
densities/currents. Similar situation also happens in ex¬ 
ternal magnetic field where the interplay between GME 
and GSE leads to a gapless collective mode called “Chi¬ 
ral Magnetic Wave” (CMW) In this Letter, we 

show that the vortical effects also induce a new wave 
mode for vector and axial density fluctuations which we 
call a “Chiral Vortical Wave” (CVW). We will derive this 
new wave equation and determine the GVW propaga¬ 
tion speed in both the hydrodynamic and kinetic theory 
frameworks. We then show that the CVW can induce a 
fermion charge quadrupole distribution from initial vec¬ 
tor density. Finally we will make predictions for possible 
implications of CVW in heavy ion collisions. 

While we will use rotating QGP as a concrete example, 
the theoretical contents of CVW is in a general fashion 
and the proposed phenomenon is universal, being able 
to be realized in various systems with effective chiral 
fermions. We note that various other chiral effects have 
been explored in systems from compact stars [55H55] . to 
Weyl and Dirac semimetals and to spin-orbit 

coupled cold Fermi gases [iiHiS] (as reviewed in mm)- It 
would be feasible and of great interest to realize CVW as 
a new and independent way of manifesting chiral effects 
in those systems. 

The Chiral Vortical Wave. — Let us start by rewrit¬ 
ing the vortical effects Q and ([^ in terms of currents 
Jl/r — 2i’^v T Ja)- 


Ja = 




U) 


( 2 ) 


Jl/r — T 




(3) 




2 


where ^j.l/r = M T M 5 - Intuitively the above vortical 
effects can be understood as follows. In the presence of 
global rotation, the underlying fermions experience an 
effective interaction of the form ~ —lj ■ S in their local 
rest frame, with S the spin of fermions. This causes a 
charge-blind spin polarization effect (as indeed found in 
other context Ea EH), namely the fermions will have 
their spins preferably aligned with us. As a result the 
right- or left-handed fermions will preferably have their 
momenta parallel or anti-parallel to u), giving rise to the 
currents in ([^. 

Let us then consider the small fluctuations of left- and 
right-handed densities on top of a uniform equilibrium 
background. For simplicity we consider oj to be con¬ 
stant and neglect fluctuations of temperature (which is 
controlled by linearized hydro equations for energy mo¬ 
mentum tensor). By combining the continuity equations 
dt-riL/R + V • Jl/r = 0 (here ul/r = ^{Jy T Ja)) with 
Eqs. ([^, one obtains, 

dtriL/R = ±^ujd^{nl/R) = ( 4 ) 

where we have set vorticity along x-direction a; = ux 
with X = a;/|a;|. Clearly there are two modes, one for 
right-handed density and the other for left-handed den¬ 
sity, that propagate in opposite direction. For later con¬ 
venience we introduce the susceptibilities for the corre- 
spending densities: Xr- = dnL/RjduLjR. 

At this point, two possibilities may occur. For simplic¬ 
ity we focus on the right-handed mode below. 

(1) The background fluid is charge neutral, i.e. /Tq = 0. 
In this case the density fluctuation Sn = Xo^h- is 
governed by the nonlinear wave equation: 

dt{Sn) + ^ ^ dxiSn^) = 0 . (5) 

47r Xq 

This takes the form of inviscid Burgers’ equation [42] 
whose (implicit) solution can be formally written as: 

„ . . ^ ujtSn\ 

Sn{x, t) = F^i^x- j ’ ( 6 ) 

with Fi{x) = Sn{x,t = 0) given by the initial density 
fluctuation. 

(2) The background has nonzero charge density /rg yf 0. 
In this case one can linearize the equation for the evo¬ 
lution of small density fluctuations on top of the back¬ 
ground density and obtain a linear wave equation: 

dt{Sn) d^{5n) = 0. (7) 

2'^ Xno 

This is just a usual wave equation describing a propagat¬ 
ing mode with a gapless dispersion relation: 

y = Va\k\ , = ( 8 ) 

271" Xr-o 


where v is the wave frequency and k = kx is the wave 
vector. This is the Chiral Vortical Wave (CVW) with 
the wave speed Vh defined above. More precisely this is 
the right-handed wave mode that propagates along the 
u} direction. The left-handed wave mode propagates in 
opposite direction to ui, with a speed given by a similar 
formula albeit replacing /ig and Xrq with the left-handed 
quantities. 

In short, the CVW found above is essentially a hydro- 
dynamic density wave arising from slowly varying vector 
and axial density fluctuations that are coupled together 
through vortical effects. Possible diffusion effects can be 
also included by adding to the LHS of Q terms like 
— DTdR){6n) where Dr and Dt are the longi¬ 

tudinal and transverse diffusion constants. If the initial 
condition of density fluctuations is very “lumpy” then 
such diffusion effects must be taken into account. 

CVW from Chiral Kinetic Theory. — Recently, the 
physics of chiral anomaly has been incorporated into ki¬ 
netic theory framework and the anomalous transport ef¬ 
fects such as the CME and CVE |47H5^ as well as the 
CMW [51| were understood in a transparent way in such 
chiral kinetic theory. It is therefore desirable to under¬ 
stand how the newly found CVW may arise in the chiral 
kinetic theory framework. 

Let us consider a rotating system of noninteracting 
right-handed (denoted by “-I-”) Weyl fermions as well as 
their left-handed (denoted by “-”) anti-particles. (The 
discussion for a system of left-handed fermions with their 
right-handed anti-fermions will be similar.) Taking a sim¬ 
ilar approach as in uni EH, we start from the equation 
of motions for these fermions in their local rest frame 

x=-±- , p = 2pxuj, (9) 

P P 

where v'G± = 1 ± p • ui/p^ and w = |a;| is the global ro¬ 
tational angular speed of the system. The corresponding 
kinetic equations can be written as: 

dtf± + X ■ d^f± p ■ dpf± = C'±[/+, /_] . (10) 

Integrating these equations and using Eq. § , we obtain: 

dt f V^/±+ f (^±'^)-9x/± 

Jp Jp P P 

+ f 2(p X u;) • dpf± = f . (II) 

J p J p 

The last term of L.H.S is zero after integration by part. 

We now examine small fluctuations in the net (vector) 
density on top of certain background equilibrium distri¬ 
bution fo±{p',T, pq). Similar to the analysis in [SH, we 
parameterize the density fluctuations as: 

f±it,x,p) = fo±{p) Sf±{t,x,p) , (12) 

Sf± = ±[dpfo±{p)] J dvd?kC^''^-'^-^^h{v,k,p) .(13) 
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where the 6f± have been expanded in Fourier modes. 
Subjecting the above into Eq. 0 , and taking a difference 
to yield the time evolution of net density, we obtain the 
following relation in linear order of fluctuations, 

ly f [dpfo+{p) + dpfo_ (p)] h{v, fc, p) 

Jp 

= k - f - [dpfo+{p) - dpfo-ip)] h{v, k,p) , (14) 
J p P 

where we have used the facts that (a) the equilibrium 
distribution /o± is a space-time independent fixed point 
of the collision kernel, (b) J^p z(p) = 0 for any z{p = 
IpI), and (c) the collision terms from fluctuations vanish 
because of charge conservation constraint. 

Let us then examine the low frequency, long wave¬ 
length limit, ^ 0 and fc —>■ 0. The hydrodynamic zero 

mode in this limit arises from 5f± —>■ ±[9p/o±(p)]iL(t, £c) 
which implies h becoming independent of p. This allows 
one to perform integrations over p in ( |14[ ), and obtain: 

^ ^ H fc) = 0. (15) 

where Xp,o is the thermodynamic susceptibility defined 
in equilibrium, = 9n/dp\T,pQ with net charge 

density n = fp '/G{fo+ — fo- )- The constant above 
is defined by C = - J^{l/p)[dpfo+{p) - dpfo-ip)] = 
^ / dp[fo+{p) — fo-{p)], and for the Fermi-Dirac dis¬ 
tribution C = po/{2tt‘^). This allows one to immediately 
identify a hydrodynamic collective excitation that prop¬ 
agates along the vorticity direction, k || w, with the fol¬ 
lowing dispersion relation: 

ty = Vn\k\ , Vh = . (16) 

271 " Xpo 

Notably, the so-obtained CVW speed Vh agrees exactly 
with that in Eq.([^. 

CVW-induced charge transport. — We now discuss in¬ 
teresting charge transport phenomena induced by the 
CVW. To be concrete we consider different initial den¬ 
sity fluctuations on top of a background medium with 
vorticity ui = wx and CVW speed Vh given by ([^ . 

Let us first consider purely axial charge density fluc¬ 
tuation in the initial condition, Ff^{x) at t = 0 (we 
have suppressed “trivial” coordinates y^z). This can 
be cast into initial conditions for right-handed and left- 
handed density fluctuations, /2. The sub¬ 

sequent evolution via the wave equations simply yields 
= ±F/^(a; =F Vht)/2. We are interested in the 
transport of vector charge density which is an observable 
quantity. This can be obtained as follows: 

{5n)Y = [FY{x-Vnt)-FY{xFVnt)\/2 

« [-d,FY{x)]Vnt (17) 

where the second line is true for small VhC This implies 
a separation effect of vector charge along the vorticity 


direction: when initial axial charge fluctuation is posi¬ 
tive and concentrated around a: = 0 (e.g. a Guassian 
form) then CVW will transport positive/negative vector 
charges toward directions respectively, leading to a 
charge dipole moment in parallel to u) with strength |dn| 
proportional to CVW speed Vh and propagation time; 
when the initial axial fluctuation is negative, the dipole 
moment flips. 

Let us then consider purely vector charge density fluc¬ 
tuation in the initial condition, F^ {x) at t = 0. Fol¬ 
lowing similar procedure, the vector charge density from 
CVW evolution is given by: 

{5n)Y = [FY{x-Vnt) + FY{x + Vnt)]l2 

^ FY{x) + [dlFY{x)l2]{Vntf (18) 

where again the second line is true for small Vot. This 
transport process leads to a quadrupole moment of vector 
charge density along ui. Consider the initial fluctuation 
to be positive and concentrated around a: = 0 (e.g. a 
Guassian form), then [d^FY{x)/2] is positive (negative) 
at large (small) |a;|, implying concentration of positive 
charges away from x = 0 toward both directions along uj. 
The resulting quadrupole |go| is proportional to (Vht)^. 

Experimental observable in heavy ion collisions. — Our 
discussions on the CVW and its induced charge trans¬ 
port effects so far are rather general. We now consider 
Its possible experimental manifestation In a concrete sys¬ 
tem, namely, rotating quark-gluon plasma created in off- 
central heavy Ion collisions. The global rotation points 
in the out-of-plane direction. In such a QGP, CVW oc¬ 
curs for each light flavors, e.g. u, d quarks and possibly s 
quarks as well and transports flavor charges toward the 
two “tips” of the QGP fireball leading to a quadrupole 
charge distribution on the transverse plane. Here we 
make a first estimate of such effect. 


<5ncvw 



FIG. 1: CVW-induced flavor charge density profile. 

Let us first quantify the quadrupole moment Qq re¬ 
sulting from CVW for a single quark flavor. We use 
the participant density from Glauber model as an ini¬ 
tial condition for the flavor charge fluctuation and study 
the dependence of on the key parameter Vh by solv¬ 
ing the GVW equation. An illustration of the GVW- 
transported flavor charge density distribution at r = 8fm 
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(with beam energy ^/s = 200GeV, impact parameter 
b = 7fm, initial time tq = 0.6fm) is shown in Figj^ 
a quadrupole pattern is evident. The quadrupole mo¬ 
ment can be obtained by integrating the density distri¬ 
bution Qq = [/ dxdy{Snf) cos{2(j)s)]/[J dxdy{Snf)\ . We 
have computed this quantity with the results: Qq ~ 
—0.03(VhAT)^ with Ar (in fmic) the propagation time 
in QGP. The numerical coefficient is for minimum bias 
events and it varies at ^ 15% level across centrality. The 
minus sign is merely due to convention of defining az¬ 
imuthal angle (j)s with respect to the in-plane direction. 

Clearly we need a plausible estimate of Vh. Us¬ 
ing the “A Multi-Phase Transport” (AMPT) model [S3] , 
our simulations suggest an initial value of lj at about 
0.5fm~^ [S4US7] . and lattice susceptibility at initial tem¬ 
perature Tq ~ 350MeV is about X/ ~ 3fm“^ [SSI jSSj, 
with both decreasing as QGP expands. Using a back¬ 
ground density /tq of 0.1 ^ lfm“^ one gets an estimate 
< Vh >~ 10“^ ~ 10“^, leading to an induced quadrupole 
\q[i\ at ~ 10“^ level. Such estimate is very sensitive to 
/io- By going to lower beam energy or by selecting events 
with large baryon asymmetry, the background /io could 
be considerably increased thus magnifying Ig^j. 


„ 1-6 
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FIG. 2: Normalized A and A elliptic flow splitting, — 
^^]/[kn|A±], for symmetric 2-flavor(2-F) and 3-flavor(3-F) 
cases. The mid-central and the peripheral correspond to for 
15 — 30% and 60 — 92% centrality class (see [60|1. 

Flavor quadrupole implies that more baryons will be 
formed on the tips than on the equator of the fire¬ 
ball. The stronger in-plane radial flow will thus trans¬ 
late the quadrupole into baryon/anti-baryon V 2 split¬ 
ting. This mechanism is in analogy to the electric charge 
quadrupole induced by CMW (33) . Suppose at the freeze- 
out, the flavor-wise chemical potential for quarks con¬ 
tains the CVW-induced quadrupole contribution (5/t/ oc 
2 (/q cos(2(^s) (with f = u,d,s). The corresponding chem¬ 
ical potential for a given type of hadron can be deter¬ 
mined from its constituent quark content, e.g. for A 
baryon SfiA oc 2{q'^ + + 9n) cos(2(/)s). We particularly 

propose to use A baryon which is electric charge neutral 
thus unaffected by possible CMW effect. We then use 
the STAR blast-wave model [60] to compute the result¬ 
ing differential flow splitting. As it is unclear how much 
the s quark mass may reduce their chiral effects, we con¬ 
sider two extreme cases: a symmetric two-flavor (2-F) 



case qn = Qn = Qn with q^ = 0, or a symmetric three- 
flavor (3-F) case q^ = dn = qQ = qh- From Cooper-Frye 
scheme it is easy to see Av 2 = oc with 

— N^)/{N^ + N^) the A-asymmetry that is 
directly related to background density /to (in analogy to 
a similar relation in CMW case [33]). The results for nor¬ 
malized flow splitting Av 2 /[\q(i\A^ are shown in Fig| 
Note that while the curves for the two centralities ap¬ 
pear close, they have rather different normalization as 
the I^qI ~ Vq strongly depends on centrality. Note also 
that the CVW predicts a particular slope for Av 2 oc A^ 
but may not exclude a finite intercept at A^ = 0 with 
either sign. Needless to say, these are crude estimates 
and a realistic hydrodynamic modeling of CVW (that 
accounts for factors like time-dependent vorticity, sus¬ 
ceptibility and diffusion) will be done in a future work. 
Given that, our results suggest that a CVW-induced sig¬ 
nal could be detected and may give indications on chiral 
effects of strange flavor. Experimental measurements of 
A and A V 2 are feasible (see e.g. |6T|), and the predicted 
flow splitting may be measured for events binned accord¬ 
ing to their baryonic number asymmetry |62] . 

Summary. — In summary, we have found a new gapless 
collective excitation in a rotating fluid system with chiral 
fermions, named as the Chiral Vortical Wave. We derive 
the wave equation for CVW and determine its speed from 
both hydrodynamic and chiral kinetic theory. We demon¬ 
strate that the CVW can induce flavor quadrupole in ro¬ 
tating quark-gluon plasma in heavy-ion collisions which 
in turn split the elliptic flow for A baryons. Such proposal 
could be tested with future experimental data. 

As a final remark, while the proposed CVW bears cer¬ 
tain similarity to the CMW, it is a completely new phe¬ 
nomenon that provides an independent way of manifest¬ 
ing chiral anomaly. In the context of heavy ion collisions, 
the vorticity lasts significantly longer than the lifetime of 
strong magnetic field and may induce more robust sig¬ 
nal. In certain spin-orbit coupled cold Fermi gases that 
could simulate chiral anomaly effects [13] , only the vortic¬ 
ity driven effects can be easily and directly induced with 
those charge neutral atoms. With such unique merit the 
CVW has its own significance and interest for the study 
of anomalous effects in a wide range of physical systems. 
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